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The  formation  of  a  cavity  in  water  is  a  conceptual  preliminary  to  the  dissolution  of  a  nonpolar  solute.  The  process  is  of  biophysical 
importance  because  of  its  assumed  pnmitive  relation  to  hydrophobic  factors  affecting  biological  structure  and  function.  In  this  paper, 
cavity  formation  by  the  isochoric  deformation  of  a  fluid  specimen  is  investigated.  Using  both  molecular  and  macroscopic  descriptions 
of  this  process,  the  volume  occupied  by  the  fluid  is  held  constant  while  the  sphencal  specimen  is  reversibly  deformed  in  differential 
steps  until  it  surrounds  a  sphencal  cavity.  The  work  of  cavity  formation  results  from  the  integration  of  energy  changes  generated 
I  directly  from  the  forced,  differential  deformations  of  the  equilibrium  fluid  structure;  the  heat  of  cavity  formation  is  the  integrated 

result  of  energy  changes  that  accompany  the  nonforced,  differential  adjustments  in  the  distribution  of  matter  which  modulate  the 
average  fluid  structure  dunng  the  process.  This  reversible  isochoric  deformation  procedure  tor  introducing  the  nonpolar  system  into 
the  pvilar  fluid  is  compared  with  the  more  familiar  melhixl  that  uses  a  reversible  coupling  parameter.  Simple  calculations  based  on  the 
isochoric  deformation  method  provide  estimates  of  the  free  energy,  internal  energy,  heat  and  entropy  of  cavity  formation  that  agree 
with  values  from  the  scaled  particle  method.^.^^ 


I.  introduction 

One  of  the  most  interesting  problems  in  bio¬ 
physical  chemistry  today  concerns  the  origin  and 
implications  of  the  hydrophobic  effect  (1-3).  In 
descriptive  terms,  this  effect  refers  to  the  tendency 
for  nonpolar  groups  to  associate  in  water.  Most  of 
the  biologically  important  issues  were  raised  in 
Kauzmann's  paper  [4]  and  the  period  since  its 
publication  has  produced  an  increasing  variety  of 
methods  and  approaches  directed  at  various  fea¬ 
tures  of  hydrophobicity  -{5].  An  assumption  in 
-mtich-of  this  work  is  that  some  fundamental  fea- 
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tures  of  hydrophobic  processes  are  reflected  in  the 
V  low  solubility  of  nonpolar  solutes  in  water  [6-11]. 
^ll  has,  nevertheless,  remained  a  matter  of  con¬ 
troversy  how  the  relative  insolubility  of  nonpolar 
molecules  in  water  is  correlated,  for  example,  with 
free  energy  constraints  on  peptide  chain  folding  in 
proteins  or  aggregation  of  lipid  molecules  in  mem¬ 
branes  .  1.  '  ..J.,, 


Of  the  concepts  which  have  evolved  regarding 
hydrophobicity,  we  state  three  propositions  which 
seem  to  have  special  significance.  These  are  (i)  the 
large  unfavorable  free  energy  of  dissolving  a  non¬ 
polar  solute  in  water  can  be  attributed  to  a  large 
entropy  decrease  associated  with  rather  special 
structural  changes  within  water  [14-16],  (ii)  the 
unfavorable  free  energy  is  not  due  to  structural 
changes  but  is  directly  dependent  on  the  structure 
of  water  [17],  and  (iii)  it  may  be  useful  to  study 
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these  processes  in  water  from  the  general  perspec¬ 
tive  of  intermolecular  force  and  energy  as  in  other 
physical  systems  [18].  The  third  statement  is  essen¬ 
tially  a  recommendation  regarding  the  other  two. 
In  pursuing  this  recommendation  here,  we  do  not 
suggest  that  it  has  not  also  been  used  in  the 
approaches  to  the  problems  referred  to  in  (i)  and 
(ii).  Further,  we  do  not  attempt  to  resolve  the 
points  of  view  expressed  in  (i)  and  (ii);  these  have 
been  thoroughly  investigated  and  defended  by  their 
proponents.  Instead,  our  purpose  is  to  study  the 
energetic  features  of  a  particular  process  that  is 
important  in  hydrophobicity.  The  process  is  that 
of  creating  a  spherical  cavity  within  a  condensed 
phase.  Cavity  formation  is  not  the  only  compo¬ 
nent  of  dissolving  a  nonpolar  solute  in  water,  but 
it  is  a  highly  significant  part  in  which  important 
features  of  water  as  a  solvent  are  found  and  is 
therefore  worthy  of  special  attention  [19-22]. 
Several  computational  approaches  to  the  thermo¬ 
dynamic  quantities  in  this  process  have  been  given; 
we  mention  here  (a)  the  use  of  scaled  particle 
theory  [19,20,23-25],  (b)  the  introduction  of  the 
cavity  by  reversible  coupling  of  an  appropriate 
hard  sphere  [26,27]  and  (c)  surface  free  energy 
relations  [28-30],  In  the  present  paper  we  will  use, 
and  assume  as  well  established,  these  earlier  theo¬ 
retical  estimates  of  positive  free  energy  and  nega¬ 
tive  entropy  of  cavity  formation  in  water;  it  is  not 
our  aim  here  to  present  alternative  calculations. 
Our  objective  is  instead  to  develop  a  description 
of  the  energetics  of  cavity  formation  through  which 
the  physical  quantities  involved  in  statements  (i) 
and  (ii)  above  may  be  explored.  This  is  done  in 
two  stages.  First,  to  introduce  in  preliminary  form 
some  concepts  relevant  to  reversible  processes,  we 
discuss  briefly  the  procedure  of  placing  a  nonpolar 
solute  in  water.  Next,  employing  these  general 
concepts,  the  creation  of  a  cavity  by  the  procedure 
of  isochoric  deformation  is  investigated,  using  both 
molecular  and  macroscopic  approaches. 

The  approaches  to  solvent  behavior  used  here 
are  stated  in  fairly  general  terms;  the  specific- 
systems  of  biophysical  interest  are  of  course  aque¬ 
ous.  As  is  well  known,  the  physical  chemistry  of 
these  systems  is  very  complex,  and  we  offer  here 
1  only  a  brief  summary  [31].  The  properties  of  liquid 

water  and  aqueous  systems  containing  nonpolar 


groups  are  primarily  determined  by  the  strong, 
highly  directional  hydrogen  bonds  among  water 
molecules.  Each  molecule  interacts  through  hydro¬ 
gen  bonding  with  up  to  four  other  molecules  that 
tend  to  be  tetrahedrally  distributed.  A  dominant 
feature  is  for  the  water  molecules  to  seek  the 
maximum  number  of  such  bonds,  even  if  some 
bonds  are  not  perfectly  formed,  i.e.,  ‘strained.’ 
The  tendency  to  maximize  hydrogen  bonding 
places  enhanced  motional  constraints  on  mole¬ 
cules  in  solvent  layers  directly  adjacent  to  a  non¬ 
polar  solute  and  is  related  to  the  strongly  negative 
entropy  of  solution.  The  energy  of  such  a  system 
is  not  given  rigorously  by  summing  interaction 
energies  over  pairs  of  molecules;  i.e.,  the  hydrogen 
bonding  interactions  are  cooperative  [32].  Thus, 
the  assumption  of  pairwise  additivity  of  inter¬ 
molecular  energies,  which  is  made  in  most  discus¬ 
sions  (including  the  present  one),  can  be  only  an 
approximation.  We  mention,  however,  that  meth¬ 
ods  have  been  given  for  obtaining  ‘effective’  pair 
potentials  that  incorporate  cooperative  many-body 
effects  [33.34]. 


2.  The  reversible  coupling  method 

As  a  preliminary  to  our  later  discussion,  we 
examine  here  a  basic  and  familiar  hydrophobic 
process,  the  introduction  of  a  small  nonpolar  so¬ 
lute,  such  as  argon,  into  liquid  water.  The  method 
developed  here  is  used  because  it  introduces  and 
illustrates  certain  points  which  are  important  for 
the  later  presentation.  The  issues  raised  in  state¬ 
ments  (i)  and  (ii)  above  are  fundamentally  of  an 
energetic  nature.  We  therefore  focus  on  the  energy 
terms  generated  when  a  single  nonpolar  solute  is 
placed  at  a  fixed,  central  location  in  a  macro¬ 
scopic  system  of  N  water  molecules  in  constant 
volume  V  at  constant  temperature  T.  In  this  sec¬ 
tion,  subscript  k  denotes  the  solute  molecule  fixed 
at  the  coordinate  origin;  subscript  j  refers  to  any 
solvent  molecule  with  position  and  Euler  orien¬ 
tation  The  configuration 

combines  the  spatial  and  orientational  coordinates 
[35].  With  k  fixed  at  the  origin,  the  intermolecular 
potential  energy  shared  between  k  and  any  j 
depiends  on  Aj  and  is  expressed  as  );  the 
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potential  between  two  solvent  molecules,  j  and 
j\  is  dependent  on  their  configurations  and  is 
given  by  Xy).  The  contributions  to  the 

internal  energy  of  the  system  made  by  intramolec¬ 
ular  vibrations,  and  by  translational  and  rota¬ 
tional  kinetic  energy,  are  considered  to  remain 
constant  throughout  the  constant-temperature 
procedure.  For  simplicity,  we  assume  pairwise  ad¬ 
ditivity  of  intermolecular  potentials.  The  method 
given  by  Kirkwood  (36),  employing  a  coupling 
parameter  k,  is  used  to  introduce  the  solute.  The 
interaction  betw'een  k  and  all  molecules  j  is  in¬ 
creased  in  reversible,  differential  steps  die,  be¬ 
tween  K  =  0  and  K  =  1 .  At  the  coupling  level  k  the 
intermolecular  potential  between  k  and  any  j  is 
KUi,^(X^)\  the  potential  between  any  two  solvent 
molecules,  Ujy(Xj,Xy),  does  not  vary  with  k. 
However,  the  potentials  of  mean  force  w^^(Xj,k) 
and  Wjy{X^,Xy,K),  which  determine  the  prob¬ 
abilities  of  occurrence  of  specific  configurations  in 
the  system,  are  functions  of  both  ic  and  configura¬ 
tion.  Thus,  the  energy  of  the  system  for  a  specified 
IC.  exclusive  of  the  vibrational  and  kinetic  energies 
discussed  above,  with  =  l/kT(kT:  Boltzmann’s 
constant  times  temperature),  is  given  by  [37] 


A* 


2(87rU’) 


/  dX.dX^ 


(1) 


With  obvious  notation,  eq.  1  may  be  expressed  as 

E(K)  =  E,,(K)  +  E^y{K)  (2) 

Formally,  each  term  on  the  right-hand  side  of  eq, 
2  has  the  form  (38) 

^  =  (3) 

I 

where  P,  describes  the  probability  distribution  in 
the  system,  over  the  manifold  of  all  possible  pair- 
configuration  energies  £|.  TTius,  £,  corresponds  to 
in  the  first  term  of  eq.  1  and  to  u^y  in  the 
second  term.  Then,  since  both  terms  on  the  right- 
hand  side  of  eq.  2  are  functions  of  the  coupling 
level,  they  may  each  separately  be  expressed  in 


differential  form  according  to 

<iE=  3£,/3ic  j  dx  0P,/3xj  dx  (4) 

As  discussed  elsewhere  [38],  the  first  term  on  the 
right-hand  side  of  eq.  4  corresponds  formally  to 
work;  infinitesimal  change  in  energies  with  con¬ 
stant  probabilities  produces  differential  work 
(-dW').  The  second  term  on  the  right-hand  side 
of  eq.  4  corresponds  formally  to  heat;  infinitesi¬ 
mal  change  in  probabilities  with  constant  energies 
produces  differential  heat  (d0).  We  apply  eq.  4  to 
each  of  the  terms  on  the  right-hand  side  of  eq.  1 
and  note  that 


,  3u^^. 


=  0 


Collecting  work  and  heat  terms,  we  have 

,V 


die 


(5) 


(6) 


/  hJ  A  \ 


I 


2{Sit^v 


2.,  .2 


J 


dXy  dX/  dx 


(7) 


Integration  of  eq.  6  over  x  yields  the  usual  defini¬ 
tion  for  the  work  of  coupling  [36],  This  approach 
to  eqs.  6  and  7  has  been  used  in  order  to  bring  out 
certain  general  points  regarding  the  coupling  pro¬ 
cedure.  First,  we  see  that  the  work  in  eq.  6  derives 
solely  from  a  summation  of  properly  weighted 
solute-solvent  interactions,  and  that  the  heat  in  eq. 
7  is  a  combination  of  contributions  from  both 
solute-solvent  and  solvent-solvent  interactions. 
Second,  we  note  that  this  coupling  procedure  is 
one  example  of  a  hypothetical,  reversible  mecha¬ 
nism  that  can  use  an  infinitesimal  change  in  free 
energy  of  a  system  to  perform  a  numerically  equal 
amount  of  work  on  the  surroundings  of  that  sys¬ 
tem.  Finally,  we  point  out  that  for  a  step  dx  to  be 
reversible,  both  energy  transactions,  the  work  in 
eq.  6  and  the  heat  in  eq.  7,  are  required;  i.e.,  each 
isothermal,  reversible  step  dx  generates  two  kinds 
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of  energy  changes  that  are  conceptually  separable 
and  have  formally  distinguishable  origins.  All  of 
these  general  considerations  appear  in  another 
method,  described  in  the  following  sections,  for 
introducing  the  nonpolar  system.  Some  of  these 
points  can  be  simplified  by  division  of  the  interac¬ 
tion  between  the  nonpolar  solute  and  polar  solvent 
into  repulsive  (hard  sphere)  and  attractive  (soft 
central)  parts;  the  simplification  occurs  when  the 
attractive  part  is  omitted.  Indeed,  it  has  been 
shown  that  the  coupling  of  a  hard  sphere  is  equiv¬ 
alent  to  the  formation  of  a  cavity  of  an  ap¬ 
propriate  size  in  the  solvent  [26,27].  Because  of  its 
conceptual  simplicity,  the  process  of  cavity  forma¬ 
tion  has  been  considered  to  be  a  primitive  model 
for  hydrophobic  phenomena  [20,27].  We  describe 
in  the  following  sections  an  additional  method  for 
introducing  the  cavity. 

3.  Cavity  formation  by  isochoric  deformation 

Our  main  objective  in  this  paper  is  to  examine 
the  energy  changes  associated  with  cavity  forma¬ 
tion  from  the  perspective  of  the  Gibbs  equation, 
which  is  the  general  expression  for  the  differential 
change  in  internal  energy  in  any  thermodynamic 
system  [39,40). 

d£  =  -dW^-p  dt> /"dS  + /I  dn -(- 4*  de  (8) 

The  notation  here  is  standard  and,  in  the  present 
and  subsequent  sections,  the  terms  on  the  right- 
hand  side  are,  in  order,  (1)  non-pv  mechanical 
work,  (2)  pv  mechanical  work,  (3)  reversible  heat 
(4)  chemical  work  and  (5)  electrical  work.  The 
systems  of  interest  here  are  closed,  electrically 
neutral  and  have  constant  volume,  so  terms  2,  4 
and  5  do  not  contribute. 

d£= -dlf'+ TdS  (9) 

Depending  on  the  context,  -dlV  and  T dS  can 
assume  various  forms.  In  the  present  and  follow¬ 
ing  sections  we  are  particularly  interested  in  the 
case  where  —dW  has  its  primitive  physical  mean¬ 
ing  and,  for  a  total  system,  is  given  by  summation 
of  the  force-displacement  terms  from  each  part  of 
the  system,  i.e.,  where  the  external  force  /  dis¬ 
places  the  internal  part  /  through  the  vector  dis¬ 


tance  ds,,  and  [39,40] 

-dW^'^l-ds,  (10) 

t 

Preliminarily  to  describing  the  procedure  in  detail, 
we  need  to  mention  some  simple  points  regarding 
cavity  formation  in  general.  First,  it  is  an  experi¬ 
mental  observation  that  stable  cavities,  even  of 
molecular  size,  do  not  form  spontaneously  in 
water.  Thus,  work  (free  energy  increase)  must  be 
done  on  a  macroscopic  specimen  of  water  in  order 
to  create  a  cavity  of  some  radius  Next  we  note 
that  a  cavity  of  radius  X,  having  been  formed  by 
intervention  of  an  external,  work-performing 
agency,  has  no  tendency  to  increase  spontaneously 
its  radius.  Thus,  in  summary,  the  free  energy 
required  to  introduce  a  cavity  is  an  increasing 
function  of  its  radius  X.  Therefore,  since  there  is  a 
functional  quantity  IF(X),  there  are  also,  by  ther¬ 
modynamic  manipulation,  the  functions  £(X)  and 
0(X)  for  energy  and  heat  of  cavity  formation, 
respectively. 

We  consider  a  spherical  fluid  specimen,  whose 
volume  V  has  a  fixed  value,  with  constant  mass  M 
and  temperature  T.  Any  point  within  v  is  associ¬ 
ated  with  a  radius  vector  r.  with  origin  at  the 
center  of  the  system.  In  our  molecular  description 
(section  4.1),  the  mass  M  is  divided  among  N 
identical  fluid  molecules  with  explicit  translational 
and  rotational  degrees  of  freedom  and  mutual 
distance-orientation  dependent  interactions.  In  our 
macroscopic  description  (section  4.2).  molecular 
properties  are  not  specified  and  the  system  is 
assumed  to  consist  of  a  continuous  spatial  distri¬ 
bution  of  matter.  A  spherical  cavity  is  introduced 
into  the  center  of  the  specimen  by  a  reversible 
process  done  in  a  particular  way.  Initially,  the 
volume  V  has  the  outer  radius  Zq  =  (3i’/47r)'''’ 
and  surrounds  a  central  spherical  cavity  of  radius 
X  =  0.  The  radius  of  the  cavity  is  changed,  by 
reversible  increments  dX,  from  X  =  0  to  X  =  t 
while  the  outer  radius  r^iX)  varies  from  (3r/47r)'^’ 
to  (t’  -I-  3v/4ir)'^'.  The  full  process  is  carried  out 
in  differential  steps  in  a  manner  such  that,  at  any 
intermediate  stage,  the  constant  volume  c  is  con¬ 
fined  between  spherical  surfaces  with  radii  X  and 
(X^ -t- 3ii/4w)'^’.  During  the  process,  only  the 
matter  and  space  between  the  spherical  surfaces 
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constitute  the  system  proper.  That  is,  for  any 
value  of  X,  the  extensive  thermodynamic  proper¬ 
ties  of  the  fluid  specimen  are  obtained  by  integra¬ 
tion  only  within  the  volume  u  contained  between 
the  spherical  surfaces  defined  by  X  and  X ).  The 
surfaces  at  X  and  ro(X)  are  variable  external 
constraints  on  the  closed,  constant  T,  v  system 
and  the  process  is  thus  the  reversible  deformation 
of  the  system  from  the  initial  to  the  final  con¬ 
straints.  Although  for  any  X  the  volume  occupied 
by  the  cavity  is  not  part  of  the  system,  forces 
between  fluid  particles  must  be  considered  as  being 
transmitted  undisturbed  across  this  volume.  This 
isochoric  deformation  requires  for  the  increase 
dX,  when  the  cavity  radius  is  X,  that  the  fluid 
matter  located  at  any  point  with  radius  vector  r  be 
displaced  a  vector  distance  df,  where 

d/  =  X-dXi7r^  (11) 

and  i  is  a  unit  vector,  located  at  r.  that  is  directed 
outward  along  r. 

Forces  originating  external  to  the  system  are 
required  to  act  at  all  points  within  the  system  in 
order  to  produce  the  prescribed  deformations.  It  is 
not  necessary  for  the  present  purpose  to  specify 
the  properties  of  this  force  field  imposed  from  the 
exterior,  other  than  to  describe  its  effects  and  the 
strength,  duration  and  direction  of  application. 
The  field  acts  simultaneously  at  all  points  in  the 
system  along  radial  directions  and,  at  any  point, 
has  only  the  strength  and  acts  only  over  the  dis¬ 
tance  required  to  produce,  in  a  reversible  manner, 
the  above-prescribed  differential  deformation  df. 
Finally,  the  external  force  field  acts  on  the  system 
only  during  the  actual  physical  displacements  in¬ 
volved  in  each  deformation  step.  That  is,  with  the 
constant  T,  u  system  at  equilibrium  surrounding  a 
cavity  of  some  radius  X,  the  force  field  is  applied 
in  the  prescribed  fashion  to  produce  the  differen¬ 
tial  displacement  df  required  at  each  point;  as  a 
result  of  this  differential  deformation  of  the  entire 
specimen,  the  external  constraints  are  moved  from 
X  and  /■g(X)  to  X  dX  and  /b(X)  +  X^  dX/ro^(X). 
The  application  of  the  field  is  then  stopped  and, 
with  the  external  constraints  held  stationary  in  the 
minutely  displaced  condition,  the  interior  of  the 
system  (still  at  constant  T)  adjusts  passively  to 
maintain  equilibrium  with  the  new  values  of  the 


external  constraints.  These  procedures  are  then 
repeated  for  each  succeeding  step  dX;  the  two 
parts  of  each  step  are  needed  in  order  that  the  full 
process,  /  dX,  be  thermodynamically  reversible. 
For  a  step  dX,  the  reversible  work  —  dlF  is  gener¬ 
ated  in  the  first  part  (deformation)  and  the  re¬ 
versible  heat  d^  arises  in  the  next  part  (adjust¬ 
ment).  The  two  parts  together  yield  the  internal 
energy  change  dF"  in  the  constant  T,  v  system  for 
the  reversible  step  dX.  We  note  here  as  an  essen¬ 
tial  point  of  comparison  between  the  coupling 
parameter  method  in  section  2  and  the  isochoric 
deformation  in  this  section  that  both  isothermal 
reversible  processes  generate  distinguishable  work 
and  heat  parts.  As  discussed  earlier,  because  the 
free  energy  of  cavity  formation  is  positive,  the 
forces  operating  within  the  equilibrium  system 
will,  on  average,  act  to  oppose  the  spontaneous 
appearance  of  a  stable  cavity.  The  function  of  the 
external  force  field  described  above  is  to  provide 
an  external  agency  which,  operating  only  during 
the  deformation  part  of  each  step  dX.  will  reversi¬ 
bly  override  these  average  internal  forces  and  re¬ 
sult  in  the  creation  of  a  cavity. 

4.  Molecular  and  macroscopic  features  of  cavity 
formation 

In  this  section,  we  use  two  alternative  ap¬ 
proaches  to  describe  the  energetic  features  of  cav¬ 
ity  formation  by  isochoric  deformation,  (1)  a  stat¬ 
istical-molecular  method  (section  4.1)  and  (2)  a 
macroscopic-phenomenological  method  (section 
4.2). 

4.J.  Molecular  aspects 

In  this  subsection  we  are  interested  generally  in 
polar  liquids  with  definable  molecular  structures 
and  interactions,  particularly  water.  The  system 
consists  of  N  identical,  fully  coupled  molecules  in 
the  spherical  volume  u  at  temperature  T.  The 
molecules  are  assumed  to  interact  in  pairwise  ad¬ 
ditive  fashion  through  hard  exchange  repulsion, 
relatively  soft  dispersion  attraction  and  powerful 
orientation-dependent  forces.  They  have  effective 
radius  r„.  For  descriptive  purposes  we  use  a 
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Cartesian  coordinate  frame  with  origin  at  the 
center  of  v.  For  the  fluid  to  surround  a 
molecular-sized  cavity  of  final  radius  \  =  t, 
molecular  centers  must  be  excluded  from  a  central 
spherical  region  with  radius  r  =  2r^  [27],  so  in  the 
initial  state  all  molecular  centers  are  confined 
within  Tp  =  and  in  the  final  state  the 

molecular  centers  are  confined  between  \  =  lr^ 
and  /•()  =  {Sr^  +  3i;/4w)'^\ 

The  configuration  of  a  molecule  /  is  denoted  by 
X,  and  consists  of  the  radius  vector  /;  from  the 
origin  to  the  molecular  center  and  the  Euler  angles 
0,.  <(>,,  'f',  defining  the  orientation  of  /  in  space 
[35].  The  Euler  angles  )  are  measured 

relative  to  local  orthogonal  axes  which,  at  any 
spatial  location  for  the  center  of  i,  have  X.  T,  Z 
axes  parallel,  respectively,  to  X.  Y,  Z  axes  fixed  at 
the  origin.  For  convenience  of  notation  in  the 
present  section,  we  denote  a  representative  pair  of 
solvent  molecules  by  /  and  j  (rather  than  j  and 
j'  as  in  section  2).  The  intermolecular  potential 
energy  shared  between  /  and  j.  «,,( A",.  A", ).  is  a 
function  of  the  two-molecule  configuration  A',,  A",. 
The  probability  of  occurrence  of  the  pair  config¬ 
uration  Af,,  Xj  in  the  system  with  cavity  radius  A  is 
a  function  of  the  potential  of  mean  force  between 
I  and  j.  which  depends  on  A",,  A',  as  well  as  on  A; 
this  potential  is  denoted  by  »',,(A)  [36,37,41]. 

The  cavity  is  introduced  at  the  center  of  the 
specimen  according  to  the  procedures  given  in  the 
preceding  section.  With  the  system  at  equilibrium 
and  A  of  any  value  between  zero  and  t,  the  force 
field  is  applied  that  is  required  to  displace  in 
isochoric  fashion  the  molecular  centers  of  all  pairs 
i  and  j  from  the  state  where  they  are  confined 
between  A  and  rp(A)  to  the  state  where  they  are 
confined  between  A-l-dA  and  r„{X)  +  dX/r^ 
( A ).  The  essential  features  of  this  force  field,  which 
acts  on  the  system  from  the  exterior,  were  dis¬ 
cussed  previously.  This  differential  deformation 
requires  that  the  centers  of  any  pair  /  and  J. 
located  with  equilibrium  pair  probability  in  any 
configuration  A',,A'^  within  v,  be  simultaneously 
displaced  in  a  radially  outward  direction  by  d/ 
and  d/^. 

d/,  =  VdA/,A^  (,2) 

d/,  =  A^  dA»yr/  (13) 


where  /,  and  are  unit  vectors  pointing  radially 
outward  from  A",  and  Xj,  respectively.  During  the 
deformation  part  of  dA,  the  action  of  the  external 
force  field  thus  ensures  well-defined  functional 
relations  between  (a)  the  cavity  radius  A  and 
cavity  displacement  dA  and  (b)  the  physical  dis¬ 
placements  of  the  molecular  centers  of  any  pair  ij 
with  specified  X^.X^.  These  displacements  of 
molecular  centers  are  done  in  a  manner  that  pre¬ 
serves  the  individual  orientations  of  i  and  j  in 
Carte.sian  space;  however,  iheir  orientation  rela¬ 
tive  to  one  another  will  in  general  change.  That  is. 
differential  displacement  of  the  molecular  centers 
of  /  and  j  along  their  outward  radial  directions, 
while  maintaining  constant  orientation  of  each 
molecule  in  the  space  of  the  Cartesian  frame,  will 
cause  differentia]  changes  both  in  relative  orienta¬ 
tion  and  in  relative  separation  of  molecular  centers. 
The  reversible  work,  performed  during  this  defor¬ 
mation  part  of  the  step  dA.  is  generated  because 
these  changes  in  relative  configuration  are  done 
against  the  forces  and  torques  operating  between  / 
and  J.  We  note  then  that  the  potential  energy  of 
the  pair  ij  changes  during  the  deformation  part  of 
a  step  d A  from  u,,(r,.4o,;  r,,io^)to  -I- d/,.uj,;r, 
+  d/^,fa),)  with  d/,  and  d/^  given  by  eqs.  12  and  13. 
The  application  of  the  field  is  then  stopped  and 
the  reversible  heat  is  obtained  when,  in  the  adjust¬ 
ment  part  of  the  step  dA,  the  distribution  of 
molecules  throughout  the  specimen  adjusts  pas¬ 
sively  to  the  new  values  of  the  constraints. 

The  discussion  above  focused  on  a  single  pair  ij 
in  a  single  pair  configuration  A',,A'^,  The  proper¬ 
ties  of  the  system  of  N  molecules  are  determined 
by  N(N-  l)/2  pairs,  with  each  pair  distributed 
statistically  over  all  pair  configurations  consistent 
with  the  external  constraints.  We  develop  here  the 
statistical  expressions  for  the  work  and  heat  con¬ 
tributions  to  internal  energy  change.  As  in  section 
2.  intramolecular  vibrational  energies  and  average 
kinetic  energies  of  molecular  translation  and  rota¬ 
tion  remain  constant;  changes  in  internal  energy 
are  then  equivalent  to  changes  in  average  potential 
energy  E  in  the  system.  For  any  cavity  radius  A 
this  will  be  [37] 


dA',  dX, 


(14) 
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The  term  multiplying  the  integral  arises  from  the 
value  at  large  N  of  N(N  —  l)/2  (the  number  of 
molecular  pairs)  and  the  normalizing  factors  for 
the  differentials  dX,  and  dX,.  The  coordinate 
integrations  are  over  the  rotational  space  (j(tf,<#>,'f') 
and  Cartesian  of  the  molecules.  The  volume 
symbol  V)^  denotes  that  the  possible  positions  of 
molecular  centers  r,  while  always  within  a  volume 
of  magnitude  u,  must  be  integrated  between  the 
radial  limits  X  and  (X^  +  3f/4?7)'^^  for  a  specified 
X.  The  Boltzmann  factor  containing  the  potential 
of  mean  force  determines  the  probability  of  the 
configuration  X,,X^  with  potential  energy 

In  principle,  eq.  14  could  be  evaluated  for  any 
X.  With  regard  to  this  equation,  a  useful  property 
of  the  isochoric  deformation  is  that  any  selected 
pair  configuration  in  the  specimen  with  cavity 
radius  X  is  derived  uniquely  from  a  specific  pair 
configuration  in  the  specimen  with  cavity  radius 
zero.  That  is,  let  (r,(0),<o,;r,(0),w^ )  be  any  specific 
configuration  contributing  the  pair  energy  u,, 
(r,(0),«,;r^(0),w^)  to  the  integration  for  £(0).  As 
the  cavity  radius  increases  from  zero  to  X,  we  may 
regard  the  isochoric  deformation  as  acting  to 
transform  the  locations  of  /  and  j  from  rjO)  and 
/;(0)  to  /;(X)  and  r^(X)  where 

r(>^)  =  (X'  +  r,(0)’)'  ■’|•(r)  (15) 

M^)  =  (^'  +  r,(0)’)"’.-(r,)  (16) 

and  i(rj  and  i(r^)  are  unit  vectors  at  the  origin 
that  point  in  the  directions  of  r,(0)  and  r,(0), 
respectively.  We  thus  have  the  transformations 

(r;(0).w,;r^{0),uij  -*(/;( X ),w, ;r, { X ),wj  (17a) 

-*  u,^(r,(X),<j,;r^(X),(a,) 

(17b) 

£(0)-*f(X)  (17c) 

As  can  be  verified  from  eqs.  15  and  16,  an  im¬ 
portant  property  of  this  transformation  is  that 

dr(X)/dX  =  (X/r(X))^  (18) 

which  gives  the  displacements  in  eqs.  12  and  13. 
This  one-to-one  correspondence  between  any 


specific  contributing  to  £(0)  and  a  unique 
contributing  to  £(X)  means  (a)  all  possible  pair 
configurations  are  accounted  for  in  the  integra¬ 
tions  for  both  £(0)  and  £(X)  and  (b)  all  pair 
configuration  energies  u,j  contributing  to  any 
£(X)  are  functions  of  X.  The  latter  point  implies 
that  eq.  14  may,  similar  to  the  case  in  section  2,  be 
expressed  formally  as  a  probability-weighted  ( £,) 
sum  over  all  pair  energies  (£,);  we  note  that  the  i 
subscript,  in  this  context,  is  an  index  of  summa¬ 
tion  for  pair  energies  and  pair  probabilities  in  [38] 

£(X)=E/',(X)£,(^)  (19) 

I 

and  the  work  -dW  and  heat  dQ  may  be  obtained 
from 

d£=  8£,/ax|  dX  H-  |££,  8£,/dX  |  dX 

(20) 

In  summary,  the  deformation  procedure  may 
be  regarded  as  a  reversible  displacement  of  all 
pairs  of  molecules  from  the  initial  manifold  of  tj 
configurations  (belonging  to  the  specimen  of  cav¬ 
ity  radius  zero)  to  the  final  manifold  of  ij  con¬ 
figurations  (belonging  to  the  specimen  of  cavity 
radius  X);  reversibility  is  achieved  by  providing  a 
‘pause’  after  each  differential  displacement  of  pair 
configurations,  in  order  to  permit  a  differential 
adjustment  in  pair  configuration  probabilities. 
These  points  are  pursued  below,  where  the  work  is 
expressed  in  the  form  given  in  eq.  10. 

An  important  application  of  eq.  14  would  be  to 
derive  the  quantity  d£  =  ffXj )  -  £(X, )  when  Xj 
=  T  and  X,  =  0.  Of  more  general  interest  here  is 
the  case  where  the  state  described  by  ffXj )  is,  by 
process  of  isochoric  deformation,  only  infinitesi¬ 
mally  removed  from  the  state  described  by  £(X|). 
Then  with  X2  =  X,  +  dX. 

£(X, +£/X)-£(X,)  =  d£=  -dW^+de  (21) 

Like  the  terms  in  eq.  1,  eq.  14  is  the  summation 
over  the  probability-weighted  potential  energy 
manifold  of  a  representative  pair.  The  infinitesi¬ 
mal  deformation  procedure,  that  forces  the  speci¬ 
men  from  £(X,)  to  £(X,  +  dX),  will  generate 
differential  changes  in  both  (1)  the  possible  poten- 
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tial  energies  of  the  representative  pair  and  (2)  the 
probabilities  of  occurrence  of  specified  configura¬ 
tions  of  the  pair.  Thus, 

d£= - -f  I  (5u,,)e  dA-,  dJT, 


requirements  of  the  isochoric  deformation,  eqs.  12 
and  13. 


8r  d/,  .2 


9/-  d/ 


(27) 


4- 

2(87rhi 

/W.f’x  r*ei,t'x 

j  u,^  S{e  dx^  dX,  (22) 

and  the  expression  would  be  evaluated  at  \  =  X,. 

The  differential  quantities  S(e  and 

are  defined  as  follows; 


We  see  for  d£  in  eq.  22  that  the  first  term  is 
the  summation  of  contributions  from  changes  in 
pair  energies  under  constant  probabilities  and  gives 
the  reversible  work;  the  second  term  is  the  sum¬ 
mation  of  contributions  from  changes  in  probabil¬ 
ities  under  constant  pair  energies  and  gives  the 
reversible  heat. 

Summarizing  these  relations. 


S{e 


)  =  {a(. 


^“,,(41 


)/9X  )  dX 


(2.3) 


I  9u,,  9/-,  9i/,, 

Su„(r,.<a,;r,.<a^)  =  .^  — -  +  — — 

8u,,  3r  du  du  i 

+  +  dX 

9r,  9X  9wj  9X  | 

(24) 

As  discussed  above,  the  displacements  in  eqs.  12 
and  13  are  not  accompanied  by  changes  in  the 
individual  orientations  of  i  and  j: 


9<j,  9<»’, 

—  =0  and-^=0 


(25) 


diV  = 


I  N' 


2(8ir^r  )^ 


+/,2(V'))i^  dX 

(28) 

x-^(c  dX^  dA,|  dX  (29) 


The  gradients  of  the  intermolecular  potential.s, 
along  the  radial  directions,  give  the  magnitudes  of 
the  intermolecular  force  in  these  directions  when 
molecules  /  and  j  are  in  the  configuration 
(r,,w,;r^,io,). 


=  -/. 


and 


Sr, 


=  -f.r 


(26) 


where  /^,  is  the  magnitude  of  the  force  from 
molecule  j,  acting  at  the  center  of  /.  in  the  direc¬ 
tion  of  r,\  has  the  complementary  meaning  for 
the  force  from  /  acting  at  the  center  of  j.  An 
individual  force,  e.g.,  will  be  positive  or  nega¬ 
tive  depending  on  whether  it  is  directed  radially 
outward  or  inward,  respectively;  this  direction  will, 
of  course,  depend  on  the  specific  X,,X^.  The 
changes  in  radial  coordinates  are  given  by  the 


The  quantities  dA’^  dA",  in  these  equations  are 
simply  the  usual  differential  elements  of  config¬ 
uration  which  must  be  integrated  over  all  config¬ 
urations  consistent  with  a  specified  X;  these 
d  A^  d  X,,  which  are  essentially  indices  of  summa¬ 
tion.  should  not  be  confused  with  the  infinitesimal 
changes  in  molecular  coordinates,  eq.  27.  gener¬ 
ated  during  the  deformation  part  of  each  dX  step. 
We  see  then  that  <3£  above  can,  in  principle,  be 
obtained  by  two  methods:  (1)  eq.  14  can  be 
evaluated  for  Xj  and  X.  separately  and  the  dif¬ 
ference  calculated  or  (2)  eqs.  28  and  29  can  each 
be  integrated  over  X  between  X,  and  Xj  and  the 
results  added  to  give  4£=  -AiV  +  AQ.  In  the 
latter  case,  —AiV  and  AQ  are  each  obtained  as 
the  limit  of  the  sum  of  infinitesimal  isothermal 
contributions,  generated  while  the  isochoric  defor¬ 
mation  procedure  is  performed  reversibly. 
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4.2.  Macroscopic  aspects 

We  derive  in  this  subsection  some  energetic 
implications  of  the  isochoric  deformation  when 
the  fluid  specimen  surrounding  a  cavity  of  radius 
X  is  regarded  as  a  continuous,  macroscopic  distri¬ 
bution  of  matter.  Our  objective  is  to  present  briefly 
certain  phenomenological  points  that  are  relevant 
to  the  molecular  discussion  above.  We  are  again 
interested  in  the  fluid  of  total  mass  M.  constant 
temperature  T.  confined  in  the  constant  volume  r 
between  the  spherical  surfaces  with  radii  X  and 
(X’  +  3t'/4ir when  a  specific  value  of  X  is  of 
interest,  this  constant  volume  will  again  be  de¬ 
noted  l>;^.  A  spherical  polar  coordinate  system  is 
convenient  in  this  case  and  the  fluid  matter  ItKaied 
at  a  point  with  radius  vector  r  within  i\  then  has 
specified  {r.9,tj)).  A  volume  element  within  is 
d('^  =  r'  sin  0  dr  d6  d<J). 

We  recall  that  the  internal  energy  of  this  sys¬ 
tem,  £(X),  can  be  changed  by  application  of  an 
external  force  field  from  £(X|)  to  £(X,),  in  a 
large  number  of  successive,  reversible,  infinitesi¬ 
mal  steps.  As  before,  each  reversible  step  consists 
of  two  parts,  (1)  a  deformation  part  in  which  the 
matter  at  all  points,  in  every  volume  element  dt’x, 
is  forced  radially  outward,  resulting  in  a  displace¬ 
ment  of  the  matter  from  any  r  to  r  +  X"dX/r‘ 
and  (2)  an  adjustment  part  in  which  passive  re- 
distnbution  of  the  fluid  among  the  volume  ele¬ 
ments  will  maintain  equilibrium  against  the 
minutely  changed  external  constraints.  In  order  to 
examine  the  differential  change  d£=£'(X  +  dX) 
-£(X)  from  the  macroscopic  perspective,  we 
focus  on  the  internal  energy  associated  with  a 
single  dcx.  Due  to  the  nature  of  the  isochoric 
deformation  procedure,  an  individual  volume  ele¬ 
ment  di'x  may  be  regarded  as  maintaining  its 
magnitude  throughout  both  parts  of  a  differential 
step  dX.  That  is,  considering  any  particular  dn^, 
the  geometric  effect  of  the  deformation  part  is  to 
compress  the  element  minutely  along  the  radial 
direction  and  to  expand  it  minutely  around  the 
tangential  directions,  while  the  matter  in  du^  and 
the  magnitude  of  du^  are  each  conserved;  thus, 
the  local  density  (mass  per  unit  volume)  in  dc^  as 
well  as  the  magnitude  of  du^  remain  constant 
during  the  deformation  part  of  dX.  However,  dur¬ 


ing  the  adjustment  part  of  dX,  while  the  redistri¬ 
bution  of  fluid  causes  no  change  in  the  magnitude 
of  any  du^.  the  local  density  within  any  dt>x  tnay 
change.  We  thus  find  as  an  essential  phenomeno¬ 
logical  point  relevant  to  both  molecular  and  mac¬ 
roscopic  methods  that,  not  only  is  the  overall 
procedure  isochoric,  it  is  also  locally  isopycnic 
during  the  deformation  part  of  any  dX.  The  inter¬ 
nal  energy  of  the  fluid  [42]  within  any  dt^  sur¬ 
rounding  a  location  r  is  u^{\)D,(\)  di\,  where 
t/,(X)  and  £',(X)  are,  respectively,  the  internal 
energy  per  unit  of  mass  and  the  local  density 
(mass  per  unit  volume)  at  the  location  r  when  the 
cavity  radius  is  X. 

The  internal  energy  of  the  specimen  is  then  [42] 

£(X)  =  /''V(X)D,(X)dix  (30) 

The  reversible  isochoric  change  from  X  to  X  +  dX 
produces  the  energy  change  d£=£(X  +  dX)- 
£(X). 

d£  =  j“‘«u,(X)£),(X)  dcx 

+  /'V(A)5£,(X)  dcx  (31) 

5tt,(X)  =  (3M,(X)/aX)  dX  (32) 

5o,(X)  =  (aD,(X)/ax)  dx  (33) 

The  quantity  6u,(X)D,(X)  de^  's>  the  work  per¬ 
formed,  under  isopycnic  conditions,  by  the  exter¬ 
nal  field  on  the  fluid  in  de^  during  the  deforma¬ 
tion  part  of  dX;  w,(X  )5Z3,(X)  di>x  is  the  energy 
change  within  dv^  due  to  the  change  in  local 
density  during  the  adjustment  part,  and  represents 
the  heat  contribution  from  dt’^  in  the  step  dX.  We 
define 

(aM,(X)/ax) dx  =  {/,(X) d/,(X)/dX) dx  (34) 

where  /,(X)  is  the  total  force  exerted  by  the  exter¬ 
nal  field  on  all  the  matter  in  divided  by  the 
total  mass  in  d  it  is  the  external  force,  per  unit 
of  mass,  that  is  required  to  produce,  in  a  reversible 
manner,  the  isopycnic  displacement  of  the  fluid  in 
di>x  at  r,  when  the  cavity  radius  is  X.  From  eq.  11, 
d/,(X)  is  the  prescribed  displacement  of  the  fluid 
located  at  r  when  the  cavity  has  radius  X  and 

d/,(X)/dX  =  (X/r)^  (35) 
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From  the  foregoing  discussion,  the  macroscopic 
counterparts  of  eqs.  28  and  29  are 

-dH/=  //'’V,(X)(AA)'Z),(\)  dc,J  dX  (36) 
d0=  [/'V(X)  3A(A)/9AdrxJ  dX  (37) 


5.  Work  and  heat  of  cavity  formation  in  water 

We  examine  here  some  quantitative  implica¬ 
tions  of  section  4.  We  restate  the  work  and  heat 
components  in  eq.  22  in  the  following  forms. 


2(87rV)- 

(.38) 

d0= - ; — f  { -/3u,,5h,  (X  ) 

2(877^)“-' 

dx,  dX,  (39) 

Where  fiw,,  is  defined  by  eq.  24  and  25,  and 
5h',,(X)  is  given  by 

6w„(X)  =  (3w,^(X)/aX)  dX  (40) 

The  probability  distribution  factors  are  identical 
in  eqs.  38  and  39;  during  any  dX,  the  relative 
contributions  from  a  specific  configuration  (  X,.  A", ) 
to  the  work  and  heat,  respectively,  are  then  de¬ 
termined  by  the  signs  and  magnitudes  of  the  fac¬ 
tors  in  brackets.  While  we  indicate  here  these 
relations  between  work  and  heat  during  a  small 
change  in  size  of  the  cavity,  we  shall  not  attempt 
a  rigorous  integration  over  X  of  eqs.  38  and  39. 
We  will  instead  assume  greatly  simplified  ap¬ 
proximations  to  these  expressions  in  order  to  de¬ 
scribe  some  of  the  contributing  effects.  These  ap¬ 
proximations  are  as  follows: 

(1)  rather  than  a  summation  of  many  infinitesi¬ 
mal  steps  (with  each  step  having  deformation  and 
adjustments  parts),  the  cavity  of  radius  t  will  be 
formed  in  one  finite  step  consisting  of  (a)  one 


deformation  of  an  assumed  initial  structure  in 
order  to  open  the  cavity,  followed  by  (b)  one 
adjustment  of  an  assumed  final  structure  around 
the  cavity;  thus,  the  work  part  emphasizes  the 
deformation  of  the  bulk  liquid  structure  and  the 
heat  part  stresses  the  structural  adjustment  to 
accommodate  the  cavity. 

(2)  Simple  molecular  distributions  will  be  as¬ 
sumed  to  represent  essential  structural  features  of 
the  initial  (43)  and  final  [24,44]  states;  thus,  the 
relevant  initial  and  final  distribution  functions  in 
eqs.  38  and  39,  respectively,  are  not  explicitly 
given,  but  are  assumed  to  produce  pre-selected 
structures  characteristic  of  water  near  4°C  (277 
K). 

(3)  Only  the  energy  changes  among  nearest- 
neighbor  molecules  in  the  first  and  second  shells 
surrounding  the  cavity  are  included;  we  thus  as¬ 
sume  in  the  one-step  integrations  of  eqs.  38  and  39 
that  and  vanish  except  in  the  immediate 
vicinity  of  the  cavity  and  are  negligible  between 
non-nearest  neighbors. 

(4)  Only  dipolar  electrostatic  energies  are  con¬ 
sidered. 

(5)  5m-,^(X)  will  be  approximated  by  an  expres¬ 
sion  |45|  derived  from  application  of  liquid-state 
theory  to  interactions  of  polar  molecules  in  a 
dielectric  medium. 

Hydrogen  bonding  among  tetrahedrally  distrib¬ 
uted  nearest-neighbor  molecules  is  assumed  in  both 
the  initial  and  final  states,  with  ccnter-to-center 
distance  t  =  2,76  A  and  dipole  moment  m  =  2.18 
Debye;  m  consists  of  a  permanent  moment  1.88 
Debye  plus  a  16%  polarization  component  [43], 
The  dipolar  energy  u  between  two  nearest-neigh- 
bors  in  these  structures  is  [43] 


-  2m-(cos  y  ),, 


T 


(41) 


where  (cos  y),,  is  the  average  cosine  of  the  angle 
between  the  dipole  vectors  of  the  two  molecules; 
(cos  is  0.33  in  the  tetrahedral  configuration 
[43).  With  these  conventions  and  the  procedures 
given  below,  the  calculations  of  Su,^  and 
for  the  initial  and  final  structures, 
re.spectively.  are  straightforward;  we  shall  omit  the 
numerical  details  and  only  briefly  outline  our 
estimates  of  work  and  heat. 
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5.  /.  Initial  structure 

We  represent  the  initial  structure  of  the  bulk 
fluid  by  a  configuration  of  eight  water  molecules. 
Two  of  the  eight  molecules,  referred  to  as  ‘central’, 
are  hydrogen-bonded  to  one  another,  and  each  of 
these  central  molecules  is  hydrogen-bonded  to 
three  ‘outer’  molecules.  Geometrically,  the  molec¬ 
ular  centers  in  this  structure  correspond  roughly 
to  atomic  centers  in  an  ethane  molecule.  The 
midpoint  of  the  bond  between  the  central  mole¬ 
cules  is  chosen  as  the  site  of  origin  for  the  cavity; 
the  inital  position  vectors  r(0)  of  the  eight  mole¬ 
cules  relative  to  this  origin  are  determined  by  the 
a.ssumed  geometry,  and  the  initial  electrostatic 
energy  of  this  structure  is  the  sum  of  seven  equal 
pair  energies  of  the  form  given  by  eq.  41.  In  the 
one-step  isochoric  deformation,  the  introduction 
of  the  cavity  of  radius  t  then  displaces  all  eight 
molecules  to  deformed  positions  r(T)  which  are 
computed  from  the  inital  positions  r(0),  using  eqs. 
15  and  16.  The  factor  (cos  in  eq.  41  is  already 
rotationally  averaged,  and  we  assume  in  our  sim¬ 
ple  calculation  that  this  factor,  as  well  as  m. 
remains  constant  for  each  pair  in  the  one-step 
deformation.  Equations  of  the  form  eq.  41,  but 
with  intermolecular  distances  appropriate  to  the 
displaced  configuration,  then  give  the  elect.'ostatic 
pair  energies  in  the  deformed  state.  The  inter¬ 
molecular  distances  after  displacement  are  readily 
derived  from  the  assumed  initial  geometry  and  the 
results  of  the  isochoric  deformation,  eqs.  15  and 
16,  Six  of  the  pair  energies  in  the  deformed  state, 
between  the  central  and  outer  molecules,  are  equal; 
the  pair  energy  between  the  central  molecules  is 
more  positive  because  the  displaceO  distance  is 
larger.  We  may  then  represent  the  change  in  total 
dipolai  energy  by  du,  and  express  this  as 

■iu,  =  iiu(c,c) 6ilM(c,o)  (42) 

where  4u(c,c)  is  the  change  in  electrostatic  inter¬ 
action  between  the  central  molecules  as  a  result  of 
the  deformation,  and  Au(c,o)  is  the  change  be¬ 
tween  a  central  and  an  outer  Molecule.  With  this 
expression  we  estimate  Au,  to  be  4830  cal/mol; 
this  is  thus  our  estimate  of  the  work  of  cavity 
formation  in  the  present  single-step  approxima¬ 
tion  to  the  integration  of  eq.  38.  For  order  of 


magnitude,  this  figure  may  be  compared  to  the 
scaled  particle  method  [20]  where  5200  cal/mol  is 
required  to  make  the  cavity  in  water  for  solutes 
the  size  of  argon  or  nitrogen.  In  the  present 
method,  we  find  that  40%  of  the  4.8  kcal  derives 
from  the  interaction  between  the  central  molecules 
and  60%  is  from  the  six  equivalent  central-outer 
interactions.  We  stress  as  a  concluding  point  in 
this  subsection  that  the  distribution  of  molecules 
resulting  from  the  one-step  deformation  of  the 
assumed  initial  structure  should  not  be  identified 
with  the  assumed  final  structure  described  below; 
in  principle,  however,  if  the  cavity  were  formed  in 
small  steps,  with  reversible  deformation-adjust¬ 
ment  in  each  step,  the  final  deformed  structure 
and  the  final  adjustment  structure  would  be  the 
same. 


Hydrogen-bonded  networks  of  water  mole¬ 
cules.  with  intermolecular  distance  t.  form  cages 
around  cavities  in  the  liquid  [24,44).  In -cold  water 
(r=  277  K).  a  closed,  nonplanar  chain  of  eight 
nearest  water  molecules  forms  the  cage  for  a  cav¬ 
ity  whose  center  is  approximately  a  distance  t 
from  the  center  of  each  molecule  [24,44],  In  this 
assumed  final  structure,  the  center  of  the  cavity, 
and  the  centers  of  any  two  cage  waters  which  are 
hydrogen-bonded  to  one  another,  then  form  an 
equilateral  triangle  of  side  length  r.  There  are  in 
this  cage  structure  eight  equivalent  pair  energies, 
with  each  assumed  given  by  eq.  41;  for  our  simple 
calculation,  the  assumption  will  be  that  only  these 
eight  pair  energies  contribute  significantly  to  the 
adjustment.  Physically,  this  is  equivalent  to  the 
a.ssumption  that  in  the  process  of  cavity  forma¬ 
tion,  eight  pair  interactions  are  transferred  from 
the  bulk-water  environment  to  the  cavity-cage  en¬ 
vironment.  In  addition  to  these  pair  energies  and 
the  assumed  final  structure  just  described,  we  also 
need  a  finite  Aw^^  in  order  to  approximate  the 
integration  of  eq.  39  in  one  finite  step.  For  a  polar 
liquid,  the  effect  on  the  h  of  two  adjacent  di¬ 
poles  i  and  j,  caused  by  the  electric  moments  in 
the  surrounding  dielectric  environment,  was  ex¬ 
amined  in  a  recent  paper  [45].  In  the  notation  of 
the  present  method  and  with  some  geometric  terms 


2.  Final  structure 
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expressed  numerically,  the  effect  on  the  of  two 
equal,  parallel,  end-to-end  dipoles  p,  and  p,. 
caused  by  removing  a  molecular  electric  moment 
from  a  molecule-sized  volume  Jr  of  the  medium 
directly  adjacent  to  ;  and  j,  is  given  by  [45| 


,  (3cos'60°  +  1)'  . 

Jm  ,,  =  (  - ^ - Jc 


Ibr" 


(43) 


where  Jr  =  4'J7/3  (t/2)',  p  =  m  cos  (109°/2)  and 
p-  =  nv  cos‘54.5°  =  /»r(cos  .  We  note  that 
this  equation  has  components  which  depend  on 
both  the  tetrahedral  geometry  (109°  )  of  the  water 
cage  and  the  equilateral  geometry  (60°)  of  the 
cage-cavity  structure.  With  the  final  structure  given 
above,  from  eq.  41  and  Jm  , ^  from  eq.  43.  we 
estimate  the  heat  of  cavity  formation  q  at  277  K. 

(/  =  -  (44) 


At  277  K.  the  heat  is  -401.''  cal/ mol  and  this  is 
our  single-step  approximation  to  the  integration  I'f 
eq.  39;  the  corresponding  entropy  of  cavity  forma¬ 
tion  is  -  14,5  e.u.  These  values  may  be  compared 
to  the  results  of  the  scaled  particle  methc'd  fi'r  298 
K.  where  the  figures  are  -4400  cal  mol  and 
-  14,7  e.u.  (20|.  Our  estimate  for  the  change  in 
internal  energy  (which  would  be  the  integrated 
form  of  eq.  22)  is  815  cal /mol,  which  mav  be 
compared  to  800  cal/mol  derived  from  the  scaled 
particle  method. 


5.. I  DepenJeme  of  work  on  cann  radius 

In  order  to  examine  the  dependence  of  the 
work  -JH'  on  cavity  size,  we  compare  some 
formal  and  practical  results  between  scaled  par¬ 
ticle  theory  and  the  present  methtxl.  The  funda¬ 
mental  expression  for  -  Jff'  in  the  .scaled  particle 
method  is  [20] 

-  Jfk'(T)  =  r{p6(A)AT}47rA' dA  (45) 

•'() 

where  p  is  the  particle  number  density  and  0'(A) 
the  conditional  singlet  probability  function  that  a 
particle  is  located  between  A  and  A  dA  when  the 
cavity  radius  is  A;  kT  has  the  usual  meaning.  The 
following  substitutions  may  be  made  in  eq.  .36: 
l)J\)  =  (mass/particle)pO(/-.A  ).  dc;^  = 


r'  dr  sin  6  dB  di)).  and  we  note  that  G(r.\)  is  the 
conditional  singlet  probability  function  that  a  par¬ 
ticle  will  be  found  at  r  when  the  cavity  radius  is 
A.  If  the  resulting  equation  for  fixed  A  is  in¬ 
tegrated  over  angular  variables,  then  expressed  as 
an  integral  over  A.  the  equation  is  then 

-  JW^(t) 

=  f[  r''’"pO[r.\)Hr.\)  drUirX'  d\  (46) 
•'()  '  •'x  ' 

where  /(r.A)  now  has  the  meaning  of  radial  force 
per  particle  required  to  produce  isopycnic  dis¬ 
placement  at  r  when  the  cavity  radius  is  A.  and  it 
will  be  recalled  that  the  outer  radius  r^^{  A  )  =  ( A’  -t- 
.3(/4ir)’  '.  With  the  volume  of  the  cavity  i;  = 
4Tr/'.3A'  and  d/;  =  47rA-  dA.  eqs.  45  and  46  may  be 
transformed 

d(  -  JH  )/dr  =pG(A)AT  (47) 

d(  ■JH  )/dr  =  r"'^'p(;(r.\)/(r.\}  dr  (48) 
•'x 

The  similarity  of  eqs.  47  and  48  is  readily  ap¬ 
parent.  particularly  when  it  is  noted  that  kT  in  eq. 
47  is  the  work  done  per  particle  at  the  surface  of 
the  cavity,  and  /(/-.A)  dr  in  eq.  48  is  the  work 
done  per  particle  in  the  spherical  layer  of  thick¬ 
ness  dr  at  r.  If  it  is  suppo.sed,  as  in  the  discu.ssion 
above,  that  most  of  the  contributions  to  the  in¬ 
tegral  in  eq.  48  are  concentrated  in  the  nearest  one 
to  two  molecular  layers  around  the  cavity,  then  we 
can  use  the  approximate  calculations  to  make 
rough  comparisons  between  the  predictions  of  eqs. 
47  and  48.  It  was  found  above  that,  at  277  K.  the 
total  work  of  cavity  formation  performed  on  the 
eight  molecules  in  the  assumed  initital  structure 
was  4830  cal/mol.  or  an  average  of  600  cal/mol 
per  molecule.  For  277  K,  the  value  of  kT  is  550 
cal/mol.  As  a  further  comparison,  we  see  in  eq.  47 
from  scaled  particle  theory  that  the  integral  over 
of  pG(A)  dt)^,  which  gives  the  total  number  of 
particles  displaced  in  making  the  cavity,  is 

-  JkF(T)/A7'.  At  298  K,  the  quantity  -JIF(t)/ 
kT  is  8.7  [20],  in  rough  agreement  with  the  as¬ 
sumed  initial  and  final  structures  in  our  simple 
approximations  for  277  K. 
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6.  Summary  and  conclusions 

Our  purposes  here  are  (a)  to  compare  the  fea¬ 
tures  of  the  coupling  parameter  method  for  intro¬ 
ducing  the  nonpolar  group  (section  2)  with  those 
of  the  method  of  isochoric  deformation  (section 

4.1)  and  (b)  to  compare  the  conclusions  of  the 
molecular  (section  4.1)  and  macroscopic  (section 

4.2)  descriptions  with  one  another  and  with  the 
energetic  concerns  of  statements  (i)  and  (ii)  in 
section  1. 

For  convenience  in  the  discussion,  the  coupling 
procedure  will  be  referred  to  as  the  'ic-process,' 
and  the  deformation  procedure  as  the  ’X-proce.ss.' 
Both  reversible  procedures  are.  of  course,  strictly 
hypothetical;  k  is  an  abstract,  dimensionless 
parameter  varying  between  zero  and  unity,  while 
X,  varying  between  zero  and  t.  is  an  actual  physi¬ 
cal  dimension  of  the  system,  the  cavity  radius. 
When  ic=  1,  the  fluid  specimen  surrounds  a  fixed 
nonpolar  solute;  when  X  =  t,  the  specimen  sur¬ 
rounds  a  fixed  cavity.  Reversibility  is  achieved  in 
each  elementary  d»c  or  dX  by  the  provision  of 
distinct,  infinitesimal  work  and  heat  components. 
The  constant-temperature  surroundings  provide 
for  the  heat  components,  eq.  7  in  the  K-proce.ss 
and  eq.  29  in  the  X-process.  There  are  interesting 
contrasts  between  the  work  components  in  the 
two  cases.  In  the  x-process,  —  dH-'  in  eq.  6  results 
from  an  infinitesimal  change  in  coupling  under 
conditions  of  fixed  configuration,  summed  over  all 
configurations;  in  the  X-process,  -dlF  in  eq.  28 
results  from  an  infinitesimal  change  in  configura¬ 
tion  under  conditions  of  fixed  (full)  coupling, 
summed  over  all  configurations.  In  each  proce¬ 
dure,  the  infinitesimal  work  is  performed  on  the 
system  by  an  appropriate  external  agency  that  can 
operate  reversibly;  i.e.,  the  work  performed  on  the 
svsiem  by  the  hypothetical  external  agency  exactly 
equals  the  increase  in  free  energy  of  the  system. 

We  have  investigated,  by  parallel  development 
of  molecular  and  macroscopic  approaches,  how 
the  concepts  of  force,  energy,  structure  and  struct¬ 
ural  change  enter  the  description  of  the  work  and 
heat  of  cavity  formation  by  isochoric  deformation. 
The  unifying  principle  for  the  two  methods  is  the 
Gibbs  relation,  eq.  8.  We  note  the  similarity  of 
eqs.  14  and  30  and  the  descriptive  corre.spondence 


between  eqs.  22  and  31,  between  eqs,  28  and  36 
and  between  eqs.  29  and  37.  For  the  work  -dfF, 
we  see  that  eqs.  28  and  36  have  the  form  given  in 
eq.  10.  The  forces  specified  in  the  force-displace¬ 
ment  terms  of  eq.  28  are  internal  in  the  sense 
discussed  at  the  end  of  section  3;  the  negative  sign 
multiplying  these  molecular  force  displacements 
indicates  the  energetic  effect  of  the  external  force 
field. 

In  section  4.1  the  fluid  ‘structure’  prevailing  for 
a  value  X  is  governed  by  e  and  ‘structural 

change'  in  the  fluid  is  determined  by  the  behavior 
of  8(e  ^’'•'''^’)/3X.  The  work  of  cavity  formation, 
the  integral  over  X  of  eq.  28.  is  derived  cumula¬ 
tively  from  minute,  forced  deformations  of  a  stat¬ 
istical  structure,  with  the  work  of  each  deforma¬ 
tion  step  dependent  on  the  static,  average  struc¬ 
ture  present  at  that  step.  The  heat,  which  is  the 
integrated  form  of  eq.  29.  is  derived  cumulatively 
from  minute,  nonforced  statistical  adjustments, 
against  stationary  constraints,  as  the  isochoric  sys¬ 
tem  of  molecules  accommodates  the  expanding 
cavity.  We  see  that  the  physical  events  in  which 
potential  energy  changes  appear  as  work  or  heat, 
respectively,  are  distinguishable  and  .separable. 

In  general,  constant  p.T  systems  are  preferred 
In  thermodynamic  descriptions  becau,se  they  re¬ 
flect  the  most  common  laboratory  situations.  Nev¬ 
ertheless.  constant  i\T  conditions  have  frequently 
been  used  in  theoretical  discussions  of  condensed 
macroscopic  systems  under  moderate  pressure  be¬ 
cause  the  pressure-dependent  contributions  to 
changes  in  thermodynamic  variables  are  small  in 
such  .systems  [27,36,37.43.46).  In  .section  .‘i.  v  ■ 
di.scussed  some  formal  and  numerical  similarities 
between  the  isiK’horic  deformation  (constant  v.T) 
and  the  scaled  particle  method  (constant  p.T). 
Further,  it  has  been  shown  by  detailed  calcula¬ 
tions  in  the  scaled  particle  approach  [20]  that  the 
pressure-dependent  term  contributes  a  negligible 
quantity  to  the  free  energy  of  cavity  formation 
and  that  this  free  energy  is  determined  only  by  the 
temperature,  the  particle  diameter  and  the  number 
density.  It  has  also  been  pointed  out  [24,44]  in  the 
scaled  particle  method  that  the  volume  available 
to  the  solvent  molecular  centers  is  the  same  before 
and  after  the  introduction  of  the  cavity;  this  is  of 
course  identical  to  the  situation  in  the  present 
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isochoric  deformation  method.  In  both  cases  then, 
the  entropy  change  accompanying  cavity  forma¬ 
tion  is  determined  only  by  the  restructuring  of 
solvent  around  the  cavity  [24,44).  Thus,  if  applica¬ 
tions  of  the  present  (constant  i\T)  methods  are 
confined  to  number  densities  appropriate  to  the 
liquids  near  1  atm,  thermodynamic  conclusions 
will  be  the  same  as  in  constant  p.  T  methods. 
Finally  we  mention  that,  although  pv  work  is 
rigorously  zero  in  both  the  tc-process  of  section  2 
and  the  X-process  of  .section  4.  the  reasons  are 
different  in  the  two  cases.  In  the  »t-proce.ss.  no  pr 
work  is  done  becau.se  there  is  no  displacement  of 
external  constraints;  in  the  X-proces.s.  for  any 
increase  dX  in  cavity  radius,  the  pv  work  done  on 
the  system  at  the  inner  surface  (radiu.s  X)  is  ex¬ 
actly  cancelled  by  the  pv  work  done  by  the  system 
at  the  outer  surface  (radius  r,,(X)). 
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